Josephson current through a correlated quantum level: 
Andreev states and tt junction behavior 



E. Vecino, A. Martm-Rodero and A. Levy Yeyati 
Departamento de Fisica Teorica de la Materia Condensada C-V, Universidad Autonoma de 

Madrid, E- 28049 Madrid, Spain 
(February 1, 2008) 

Abstract 

The Josephson transport and the electronic properties of a quantum dot 
characterized by a single level coupled to superconducting leads is analyzed. 
Different approximations are used and compared: the mean field approxi- 
mation, the second-order perturbation theory in the Coulomb interaction 
and the exact diagonalization in the zero band-width limit. The system ex- 
hibits a rich behavior as a function of the relevant parameters. We discuss 
in detail the conditions for the observation of tt junction behavior and the 
effect of Coulomb interactions on the Andreev states. 

I. INTRODUCTION 

The observation of the Kondo effect in semiconducting quantum dots [1] has opened a 
new area of research in which electronic transport through a nanoscale strongly correlated 
system can be studied under controlled conditions. More recently, this effect has also 



been observed in carbon nanotubes coupled to metallic electrodes, a system that can 
behave in many respects as a quantum dot [2]. These type of devices provide an almost 
ideal system to test different theoretical predictions. In this direction, a great theoretical 
interest has arisen in connection with the possibility of replacing the normal electrodes 
by superconducting ones. In fact, carbon nanotubes connected to superconducting leads, 
that are already being produced, can provide a physical realization of such a system [3]. 

In these type of systems an issue of fundamental relevance is the interplay between 
electron correlation effects and Andreev reflection processes, which provide the basic mech- 
anism for transport between normal and superconducting regions. This interplay has been 
analyzed by several authors for the case of quantum dot coupled to a normal and a super- 
conducting electrode [4,5]. There have also been some works addressing the problem of 
the electron transport mediated by multiple Andreev reflections through a resonant level 
between two superconducting electrodes [6]. 

When both electrodes are superconducting a more basic question is how the electron- 
electron interactions in the dot would affect the Josephson current. This issue has re- 
ceived considerable attention in recent years [7-11]. In particular, the discussion has been 
centered to a large extent around the appearance of a 7r-junction behavior induced by 
electron-electron interactions. The 7r-junction behavior, which consists in a reversal of 
the sign of the supercurrent under certain conditions, was flrst pointed out by Kulik [12] 
when analyzing the Josephson tunneling in the presence of magnetic impurities, and dis- 
cussed afterwards by several authors [13-15]. This issue is also related to the interplay 
between superconductivity and magnetism. 7r-junction behavior has in fact been recently 
observed in S-F-S Josephson junctions [16]. 

The existing theoretical works analyzing the Josephson effect in a single correlated 
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level coupled to superconducting leads have adopted either a mean field description [8] 
or are restricted to some limiting situations [7,9]. In Ref. [7] the current was obtained 
to the lowest order in the tunneling coupling, thereby neglecting the important physics 
associated with the Andreev bound states. On the other hand, in Refs. [9,10] the hmit of 
infinite Coulomb repulsion was considered. 

At this point we believe that further work is needed to understand the physical behav- 
ior of this system for a broad parameter range. This behavior is actually rather complex 
due to the several energy scales involved. In particular, as the Andreev bound states 
play a crucial role in the transport properties, a detailed analysis of electron correlation 
effects on these states seems desirable. In this paper we present a theoretical study of 
the Josephson effect in a single level quantum dot coupled to superconducting electrodes. 
Different approaches are considered. We first study the zero band-width limit in which 
the problem can be exactly diagonalized. Then, as an approximation to the full model we 
use a diagrammatic expansion of the self-energy associated with the Coulomb interaction 
in which the coupling to the superconducting leads is taken into account up to infinite 
order. This last ingredient is important for a proper description of the Andreev bound 
states. 

In section II, we introduce the theoretical model and discuss the diagrammatic ap- 
proximations for the electron self-energy. In section III we present a simplified version of 
this model in which the Coulomb interaction is replaced by an effective exchange field. 
This simple model already describes the vr-junction transition and allows to understand 
the behavior of the Andreev states under a finite magnetization in the dot. Section IV is 
devoted to the analysis of the zero band-width limit in which the model becomes equiv- 
alent to a finite system which can be solved exactly. The exact solution is then used as 
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a test of the approximations used for the electron self-energy. It turns out that the zero 
band-width limit can describe rather accurately most of the properties of the full model. 
The results for the full model are presented in section V. We first discuss two opposite 
limiting situations in which the coupling of the dot to the leads is either small or large 
compared to the superconducting gap. As in the zero band-width limit, in the case of 
weak coupling we show that the problem can be solved exactly. We also present some 
numerical results for the intermediate regime. The paper is closed by a brief summary of 
the main conclusions. 



II. MODEL AND THEORETICAL APPROACH 

We describe a small quantum dot connected to superconducting electrodes by means 
of a modified single-level Anderson model, given by 

H ^ Hl + Hr + J2 ^oUa + Uh^hi + Ht, (1) 

where = Cq^Coct, Hl and Hr represent the uncoupled superconducting leads; Ht = 
^keL,R;ato,kCocr^k,a + h.c. describing the coupling between the dot level and the leads. 
Within this model the dot is represented by a single spin degenerate level with a repulsive 
Coulomb interaction described by the [/-term in Eq. (1). We shall assume that the 
superconducting leads are well described by the BCS theory with a superconducting gap 
A and that there is a fixed superconducting phase difference (f) — (pL — <f>R between both 
electrodes. 

The relevant quantities like the current and the spectral densities can be expressed in 
terms of non-equilibrium Green functions [17]. For the description of the superconducting 
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state it is useful to introduce spinor field operators [18], which in a site representation are 
defined as: 



( \ 



(2) 



Then, the difi^erent correlation functions appearing in the Keldysh formalism adopt the 
standard causal form: 



(3) 



where T is the chronological ordering operator along the Keldysh time contour. The 
labels a and (5 refer to the upper {a = +) and lower {a = — ) branches on this contour. 
The functions Gfj~, which can be associated within this formalism with the electronic 
non-equilibrium distribution functions, are given by the (2x2) matrix: 



Gtr{t,t') 



> 



(4) 



[<c]^{t')cl{t)> < Cj^{t')cl{t) 

In terms of the Fourier transform matrix elements of Gfj~{t,t') one can write the 
charge and the induced order parameter on the dot, as well as the Josephson current: 



1 f°° 



(5) 
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d(jj[G+Q (a;)]2i. 



Il,r = t / ^'^Tr (j^io.feG'Jo (a;) - a^tk,oG^^ki^) 



(6) 
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where Il,(r) denotes the current between the left (right) lead and the dot, cr^ is the usual 
Pauli matrix, and io,fe is the hopping matrix in the Nambu representation 



io,k — (ik,oj — 



to,k 

-t 



(8) 
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For the zero voltage case the calculation of the different (a;) is particularly simple 
because the following relation holds: 



(9) 



where /(a;) is the Fermi distribution function, and G""'"^ are the advanced and retarded 
Green functions. Therefore, the relevant quantity to be determined is the dot retarded 
Green function, which in a Nambu 2x2 representation adopts the form 



(10) 



where F^^/j are the tunneling rates given Fj 



T^PFio,kg''ik,o, with g'l^ 
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—uj/y/A'^ — up' and = [g. 



Ii)* 



A./^/ A?' — uPe^'^'' with k — L,R and pf is the nor- 



mal density of states at the Fermi level (the chemical potential of the superconducting 
leads is taken as zero). The self-energy E''(tf) takes into account the effect of Coulomb 
interactions. To the lowest order in U this is given by the Hartree-Fock Bogoliubov ap- 
proximation: (E'-)n <ni>, (EO22 ^-U <n^ >; (EO21 = (^0^2 = U{cl^cl^). We 
shall discuss in II A how correlation effects beyond this mean field approximation can be 
included. 

This model has been analyzed within a mean field approximation in Ref. [8]. In that 
work the mean field solution was simplified by neglecting the induced order parameter 
in the dot and imposing self-consistency only in the dot magnetization. However, the 
complete Hartree-Fock Bogoliubov solution requires the self-consistent determination of 
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both the diagonal and non-diagonal charges in the dot. It should be noticed that the self- 
consistent determination of the induced complex order parameter is in principle necessary 
in order to ensure current conservation [19]. 

A. Inclusion of correlation effects 

In order to go beyond the mean field approximation, we consider the diagrammatic 
expansion of the self-energy in terms of one-electron propagators in Nambu space. In Fig. 
1 we show the corresponding second order diagrams in the electron-electron interaction. 
Due to the appearance of the anomalous propagators G12 and G21 in the superconducting 
state, there are additional diagrams to the one contributing in the normal state, labeled 
by 11(a) in Fig. 1, corresponding to the interaction of a quasi-particle with an electron 
hole pair with opposite spin. 

The proper choice of the unperturbed one-electron Hamitonian over which the dia- 
grammatic expansion is performed is an important issue. For the normal Anderson model 
in the symmetric case (eo = —U/2) the Hartree approximation which renormalizes the 
dot level as eg + U/2 is the adequate starting point for the perturbation theory as it 
automatically warrants charge consistency between the perturbed and the unperturbed 
situation. However, in a non-symmetric case perturbation over the Hartree field yields 
pathological results close to half-filling (see [21]). As discussed in Ref. [21] a better choice 
is to define an effective dot level eg// in such a way that charge consistency between the 
effective and the interacting problems is achieved. The natural extension of this scheme 
to the superconducting case is to impose also consistency in the non-diagonal charge riu 
by introducing an effective local pairing potential Ae// in the unperturbed Hamiltonian 
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[5]. 



With these definitions, the dot Green functions in the unperturbed effective problem 
are given by (hereafter we omit the subscript 00 in the dot Green functions) 



D{u;) 



\ 



\ 



11) 



where D{u!) is the corresponding determinant. Notice that the Andreev states in the 
unperturbed problem are determined by the condition 



D{uj) = 0, 



(12) 



which, in the general case, can have up to four solutions. 

The different contributions to the self-energy represented in Fig. 1 can then be written 
in terms of the one-electron propagators in the following way 

[/2 



dei j de-2 j 
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Notice that /2m and G^^^ ^ /2m in the above equations correspond respectively 



to the occupied and unoccupied states in the dot spectral density. 

The total self-energy is obtained by adding the contributions labeled by (a) and (b) in 
the above equations. The remaining self-energy components (^^22' ^^^d ^12^'') ^'^^ obtained 
by similar expressions. In the non-magnetic case the diagonal self-energy components are 
related by E22^'*(a;) = —T!i^\—uj). This relation does not hold in the general magnetic 
case, except when eo = —U/2 for which Tl'^^\uj) — E22^^(a;). 

In order to extend the range of validity of the second order self-energy an ansatz which 
interpolates correctly between the limits [//F — > and U/T —>■ 00 can be used [5]. Within 
this interpolative approach the matrix self-energy is given by the following expression [5] 



and E(2) is the second order self-energy. 

This interpolative scheme has been used in many different contexts involving strongly 
correlated electrons [20-22]. Notice, however, that for the electron-hole symmetric case 
the interpolative self-energy reduces to the second order one. 

III. TRANSITION FROM NORMAL TO tt JUNCTION BEHAVIOR: A SIMPLE 

MODEL 

An exactly solvable model which can describe the transition from normal to tt junction 
behavior is obtained by substituting the interacting region by a single site with a local 
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where 
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ep + (1- < n >)U - eeff 
U'^ <n> (1- < 77, >) 
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exchange field E^x, in such a way that — eo + Eg^ and ei — Sq — Eg^- It should be 
noticed that this model is formally equivalent to a mean field solution of Hamiltonian (1) 
with the prescription Eg^ = U{ni — n^)/2. 

This model exhibits in general four bound states inside the superconducting gap (An- 
dreev states) which give the dominant contribution to the current. This is determined 
by the derivative of the states below the Fermi energy with respect to the phase. The 
position and spectral weight of these states can be obtained from the retarded Green 
function (Eq. 10). In the limit A << F and when F^ = F^^ = F, the Andreev states can 
be determined analytically by the expression 



It can be verified that in the limit X — > Eq. (18) recovers the well known expression 
for the Andreev states in a single channel contact of transmission r = 1/(1 + E'^), i.e. 
a;/A = ±^1 -Tsin^ 0/2 [23]. 

The qualitative behavior of the Andreev states can be easily analyzed in the particular 
case £■ = 0, in which the position of the four states is given by 



This expression clearly shows that the effect of a finite exchange field is to split the non- 
magnetic Andreev states at ±A cos 0/2 into four states. The transition to the tt junction 
behavior is associated with the progressive inter-crossing of the two "inner" bound states. 




2 



(18) 



where E = eo/2F, X = Eg^/^T and = - E"^ and 5(0) is given by 





(19) 
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This crossing is illustrated in Fig. 2 for a generic situation. For small exchange field 
(Fig. 2 a), the inner states do not cross and the system is in the state. As X increases, 
the inner states cross at tt and thus the current-phase relation exhibits a sign change 
around (p — tt. Eventually, when the crossing between the inner states is complete the 
whole current-phase relation changes sign and the system is in the tt state. The state in 
the intermediate region is conventionally designed as 0' or tt' depending on the relative 
stability of the minima of the inter-crossing states. 

The boundaries between the different regions are straightforwardly obtained from Eq. 
(18): the curves X = E, X = {E + V3 + 4E2)/3 and X = ^/l + E^ correspond to the 
— 0', 0' — tt' and n' — n boundaries respectively. The full phase diagram of this simple 
model is illustrated in Fig. 3. 

IV. ZERO BAND- WIDTH LIMIT 

A special limit in which the model given by Eq. (1) becomes exactly solvable is the 
case where the band-width of the electrodes tends to zero. This is equivalent to neglect 
the high energy excitations in the superconducting electrodes, substituting the frequency 
dependent electrodes self-energy by an effective diagonal and non-diagonal level in Nambu 
space [24]. In this limit the semi- infinite leads connected to the dot are replaced by an 
effective single site and the terms Hl,r describing the leads in Eq. (1) are given by 

HL,R = J2^L,Rci^RaCL,Ra + ^L,RCl,R^Cl,RI + A^^rC^^C^^. (20) 
a 

The system thus becomes a sort of "superconducting molecule" with a finite number 
of electron configurations which can be diagonalized exactly. This solution can be useful 
as a test of the accuracy of the proposed approximations for the infinite system. In spite 
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of its simplicity this limit also gives a qualitative description of the behavior of the full 
model, as will be shown below. 

The exact solution in this limit is obtained by considering electronic configurations 
with all possible total number of electrons, the pairing interaction connecting configura- 
tions that differ in two electrons. The eigenstates can be classified into those arising from 
even or odd number configurations. In the even case the ground state has zero total spin, 
while in the odd case the ground state is two-fold degenerate corresponding to 5'^; = ±1/2. 

In Fig. 4 we show the evolution of the ground state energy for the even and odd cases 
for fixed eo and increasing U. As can be observed, in the even case the ground state 
energy exhibits always a minimum for = 0, while in the odd case the minimum appears 
for (f> = TT. This is to be expected from the fact that the system magnetization is finite 
in the odd case. For small values of U the even case is more stable for all values of the 
phase (Fig. 4 a), while the opposite situation is found for large enough U (Fig. 4 d). 
In the intermediate region (Fig. 4 b and c) the even and odd energy levels cross. This 
level crossing corresponds to the transition between and tt junction behavior in the full 
model. The tt state thus corresponds to a situation in which the dot adquires a finite 
magnetization. 

The exact ground state energy of this simple model can now be used to check the 
accuracy of the different approximations discussed in Sect. II. For the sake of simplicity 
we shall restrict this comparison to the eo = —U/2 case. In Fig. 5 we show the comparison 
between the exact, the mean field, and the second-order self-energy results for the ground 
state energy for different values of C//A. At small values of C//A (Fig. 5a), the energy of 
the mean field approximation lies slightly above the exact result both in the non-magnetic 
and in the magnetic case. When increasing U/A (Fig. 5b), the non- magnetic mean field 
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solution increasingly deviates from the exact result. In contrast, in the magnetic case, 
the deviation between the exact and the mean field solution first increases in the small 
U range while it progressively decreases for large U. This is due to the fact that the 
exact ground state corresponds to a localized spin at the central site for [/ — > oo, which 
is correctly described by the magnetic mean field solution. 

The inclusion of the second-order self-energy substantially improves the results of the 
mean field approximation both in the magnetic and in the non-magnetic cases for small 
and moderate values of U. In the non-magnetic case the improvement is considerable in 
the whole range of values of U (Fig. 5a and 5b). On the other hand, in the magnetic 
case the improvement due to the inclusion of correlation effects progressively becomes less 
important as the mean field tends to the exact solution in the large U limit. 

Finally, it is interesting to analyze the phase- diagram of the model in the zero band- 
width limit, which turns out to contain the essential features of the full model with the 
four types of solutions: 0, 0', n' and tt appearing in the previous section. In this discrete 
model the different phases can be identified by the phase-dependence of the ground state 
energy. 

The exact phase diagram for the zero band-width limit is shown in Fig. 6a. Only 
the region U > and e < is shown as for all other regions the only possible phase is 
type. Roughly, the tt state is found for U/T and — eo/F sufficiently large, while the state 
appears either for U/T or —eo/F sufficiently small. The 0' and the tt' behavior is found 
in the intermediate regime. It is interesting to notice that for —eo/F < 1 and U/T — > oo 
the state is always the more stable. An extrapolation of this result to the full model 
would imply the absence of tt junction behavior in the so called mixed valence regime, as 
has been predicted in Ref. [9]. 
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For comparison, the mean field diagram is also shown in Fig. 5b. Notice that although 
the overall behavior is captured in the mean field solution, correlation effects displace the 
TT region to much larger values of U /T and — eo/F. Neglecting the induced order parameter 
in the dot as done in Ref. [8] does not quahtatively change the phase diagram although 
the full self-consistent solution is somewhat closer to the exact result. 

V. RESULTS FOR THE FULL MODEL 

The mean field approximation for the full model has been analyzed in Ref. [8]. Al- 
though in this work self-consistency in the induced order parameter was neglected, its 
effects on the total energy is small as already discussed for the zero band- width limit. In 
the present section we will concentrate in discussing the effects of electronic correlations 
beyond the mean field solution. We will consider the electron-hole symmetric situation 
(eo = —U/2) in which electron correlation effects are expected to be more important. 

For a fixed value of U, the physical behavior of the model is controlled by the di- 
mensionless parameter A/F. Is is interesting to analyze in detail the two opposite limits 
A/F > 1 and A/F < 1 

A/F ^ 1 hmit: This situation corresponds to a dot very weakly coupled to the leads. 
In this limit the exact solution can be obtained due to the fact that the problem can be 
mapped into a two-level Hamiltonian describing the dynamics of the Andreev states. This 
can be shown by first considering the non-interacting {U — 0) situation. In this case, the 
spectral density of the dot exhibits bound states at energies tOg = ±2Fcos0/2 and the 
continuous part becomes negligible. Consequently, the retarded Green function at the dot 
can be simply written as 
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U —Us 

— o;, cu 



(21) 



indicating that the corresponding effective Hamiltonian is given by H'^f'f = a;5[c|Cj^ + c|c|]. 
When introducing a finite U, the diagonal level at —U/2 is canceled by the Hartree 
term and the remaining part of the Hamiltonian associated with the Coulomb interaction 
becomes U{nf — 1/2) (n^ — 1/2) — t//4. This term vanishes for the case of an empty and 
a doubly occupied dot, while it yields an energy —U/2 for the single electron case. As 
a consequence, the ground state of the system is either the symmetric combination of 
empty and doubly occupied configurations (with energy -cUs) or the doubly degenerate 
single occupied state (with energy —U/2). Thus, the transition to the magnetic tt state 
takes place when U/2 > 2r. 

A/r -C 1 limit: In this opposite limit one would expect to recover gradually the prop- 
erties of a normal system. In particular, for [/ > F the features associated with the Kondo 
effect should emerge. Fig. 7 shows the dot spectral density in this regime for increas- 
ing values of U obtained in the second-order self-energy approximation. For comparison 
the spectral density in the normal state is also shown. As can be observed, the spectral 
density is similar to the one found in the normal state except for the superimposed fea- 
tures in the gap region. The overall shape evolves as in the normal case from a single 
Lorentzian broad resonance for [/ < F to the three peaked structure characteristic of 
the Kondo regime when U > T. In this regime the relevant energy scale is set by the 
Kondo temperature Tk which essentially measures the width of the Kondo resonance in 
the normal state. Within the second-order self-energy approximation Tk ~ F/(1 — ao), 
where ao — |^(0). In the symmetric case 
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which coincides with the perturbative resuh of Ref. [25] for the Anderson model. Although 
this approach fails to give the exponential decay of the Kondo temperature with U it gives 
a rather good description of the spectral density in the moderate U/T range [26]. 

The coexistence of the Kondo and Josephson effects is to be expected as far as > A. 
When U is further increased the system should evolved into the magnetic tt state with 
the suppression of the Kondo effect. 

Let us analyze the self-energy in this limit. The effective one-electron problem in this 
case is characterized by the presence of Andreev states which, from Eq. (11), are located 
at a;, = A^l -Tsin^ 0/2, where r = ^T'^/{elff + 4r2) is the normal transmission in 
the effective problem (in the electron- hole symmetric case considered here e^ff = and 
T — 1), like in the case of a single channel point-contact. The weight of the Andreev 
states at the dot site decreases as A/F according to the expression A| sin0/2|/4r. Also 
the induced order parameter tends to zero as A/F in this limit. 

The electron self-energy can then be evaluated retaining only terms of order A/F. The 
contributions labeled as 11(b) and 21(a) can be altogether neglected as they involve more 
than one anomalous propagator and are thus of order (A/F)^. On the other hand, from 
the general expressions for the second-order self-energy, Eqs. (13,16), it should be noticed 
that there are three types of contributions to E^^^ : one involving only the discrete part of 
the one-electron spectral densities, another one involving only the continuous part; and 
terms in which both the localized states and the continuous spectrum are mixed. The first 
of these contributions is of order (A/F)^ and can be neglected. The resulting expression 
for up to first order in A/F for \uj\ < A is 
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sin(0/2)| + 2(l + -)cos(0/2) A. 



a; 



(23) 



(24) 



These expressions for the self-energies allow to determine the renormalization of the 
states inside the gap due to the interactions. For moderate values of U {U < lOF) 
the renormalized states have approximately the same phase-dependence as in the non- 
interacting case (i.e. ~ cos (0/2)) but with a narrower dispersion given by uJsi^) — 
A(l - (U/Uof), where {Uo/Tf = (r/A)7rV(27r + 2). For larger values of U the phase- 
dependence of the states starts to deviate from this simple law. This is illustrated in Fig. 



The evolution of the renormalized Andreev states indicate that the critical currents 
are suppressed as ~ (1 — {U /UqY) in this limit (for moderate values of U). 

One can summarize the results for the full model in the symmetric case by discussing 
the phase diagram shown in Fig. 9. In this figure we compare the results of the mean 
field and the second-order self-energy approximations for the critical C/, [4, defining the 
transition to the tt state. As can be observed both approximations yield the same result 
in the F/A <C 1 limit. A close inspection of this regime, illustrated by the inset in Fig. 
9, shows that the transition takes place for Uc ~ 4r when F ^ in agreement with the 
exact result. As F/A increases the predictions of the two approaches start to deviate. 
The mean field predicts an almost linear relation between Uc and F although with a larger 
slope than in the F <^ A limit. On the other hand, the second order self-energy predicts a 
faster increase of Uc with F. Notice that in the normal state (A — > 0) there should be no 
transition into a magnetic state for finite U . We also show as a reference the line where 
Tk in the normal state, given by the second-order self-energy approximation, matches A. 



8. 
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It can be observed that the criterium TV ~ A for the transition to the tt state is a rather 
good for r sufficiently large. 

VI. CONCLUSIONS AND FINAL REMARKS 

We have presented a theoretical analysis of the Josephson transport through a strongly 
correlated level coupled to superconducting electrodes. The analysis has been specially 
focused in the effect of electron correlations on the subgap Andreev states. These states 
determine to a large extent the current-phase relation in this system. The transition to 
a TT state can be understood as a result of the intercrossing of the subgap states induced 
by an increasing Coulomb interaction. Within this model this transition corresponds to 
a truly quantum phase transition in which the ground state becomes degenerate having 
a locahzed magnetic moment at the dot level as already noticed in Ref. [10]. It is worth 
noticing that this situation cannot exist in the absence of superconductivity. In fact, this 
behavior can be traced to the suppression of the Kondo effect (in which the electrons 
at the dot level couple to a singlet) due to the absence of low energy excitations in the 
superconducting leads. 

In the present analysis we have used different approximation methods. In order to get 
insight on the behavior of the Andreev states in the transition to the tt state we have first 
analyzed a simple mean field model in which the electron interactions are represented by 
a local exchange field. We have also studied the zero band-width limit which allows for an 
exact diagonalization. It has been shown that the study of this limit already illustrates 
the different types of behaviors that can be found in the full model. These results have 
also been used to determine the accuracy of the self-energy approach, showing that it 



18 



improves considerably the results of the mean field approximation for moderate values of 
the Coulomb interaction. Finally, we have presented results for the full model in different 
regimes. In the limit A ^ F we have shown that the problem can be solved exactly, its 
dynamics being described by a two level Hamiltonian corresponding to the Andreev states. 
On the other hand, for F » A there is a coexistence of Kondo and Josephson effects for 
Tk > A, the main effect of electron correlations being included in a renormalization of 
the critical current. 

The present work constitutes a first step in the study of the transport properties of a 
quantum dot coupled to superconducting leads in a non-equilibrium situation, i.e. with 
an applied bias voltage. This would allow to analyze an experimental situation like the 
one of Ref . [3] . Work along this line is under progress. 
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FIGURES 
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FIG. 1. Diagrams contributing to the second-order self-energy in the superconducting state. 
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FIG. 2. Bound states within the energy gap and current-phase relation for the simple model 
discussed in Sect. III. For eo/2T = —0.5 and Eexf^I" = 0.25 (top panel), 0.75 (middle panel), 
1.50 (lower panel). 
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FIG. 3. Phase diagram for the simple model of Sect. III. 
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FIG. 4. Energy levels for the model in the zero band-width limit with ti = tn = 1.2A, 
e/A = -10 and U/A = 11 (a), 13 (b), 15 (c) and 18 (d). Dotted and fuU lines correspond to 
the even (ground and first excited states) and odd cases respectively. 
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FIG. 5. Comparison of the ground state energy in the zero band-width hmit obtained by 
different approximations: exact diagonahzation (full lines), mean field approximation (dotted 
lines) and second-order self-energy (dashed lines). The values of the parameters are e = —U/2, 
t^ = tR = 1.2A with C//A = 2.5 (a) and [//A = 10 (b). 
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FIG. 6. Phase diagram in the zero band- width limit with ti = tn = A obtained by exact 
diagonahzation (a) and within the mean field approximation (b). 
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FIG. 7. Density of states for the symmetric case in the regime T » A. The values of the 
parameters are A/F = 0.05 and U jV = 5 (a), 10 (b) and 20 (c). The dotted lines indicate the 
corresponding results for the normal state. 
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FIG. 8. Renormalized bound states within the gap in the regime F A with C//F = 0, 10, 
14, 18 and 20. 
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FIG. 9. Phase diagram in the symmetric case for the full model. The full and the dashed 
lines correspond to the onset of the tt state within the second-order self-energy and the mean 
field approximation respectively. The dash-dotted line corresponds to Tk = A within the 
second-order self-energy approximation. Inset: closer view of the F <C A region showing the 
behavior of the mean field approximation. The dotted lines indicate the different slopes in the 
r ^ A and T > A regimes. 
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